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Abstract 

We first show that the quasiclassical limit of the squared eigenfunction symmetry 
constraint of the Sato operator for the mKP hierarchy leads to a reduction of the 
Sato function for the dispersionless mKP hierarchy. The constrained dispersionless 
mKP hierarchy (cdmKPH) is obtained and it is shown that the (2+l)-dimensional 
dispersionless mKP hierarchy is decomposed to two (l+l)-dimensional hierarchies of 
hydrodynamical type. The dispersionless mKP hierarchy with self-consistent sources 
(dmKPHSCS) together with its associated conservation equations are also constructed. 
Some solutions of dmKPESCS are obtained by hodograph reduction method. 

PACS number: 02.30.IK 



1 Introduction 

Recently much attention has been focused on the study of dispersionless integrable hier- 
archies which have important applications from complex analysis to topological field theory 
(see [1-16]). In [2, 3, 6, 7, 9], a standard procedure of dispersionless limit of integrable 
dispersionfull hierarchies is proposed. In this procedure, dispersionless hierarchies arise as 
the quasiclassical limit of the original dispersionfull Lax equations performed by replacing 
operators by phase space functions, commutators by Poisson brackets and the role of Lax 
pair equations by conservation equations (or equations of Hamilton- Jacobi type ). Several 
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strategies have been proposed to investigate these hierarchies, such as the twistorial method 
[8, 9, 17], the hodograph reduction method [6, 7, 18] and the quasi-classical d-method [13- 
16]. Lately in [16], using the quasi-classical 9-method, the authors studied the symmetry 
constraint of the dispersionless KP (dKP) and 2DTL hierarchies and some important reduc- 
tions of these hierarchies are shown to be nothing but results of symmetry constraints. This 
offers us another important way to construct reductions of dispersionless hierarchies. 

The integrable equations with self-consistent sources are another type of important in- 
tegrable models in many fields of physics, such as hydrodynamics, state physics, plasma 
physics, etc (see [19-27]). In our general approach proposed recently in [22-27], the con- 
strained integrable hierarchy may be regarded as the stationary system of the corresponding 
integrable hierarchies with self-consistent sources. From this observation, the integrable hi- 
erarchies with self-consistent sources and their Lax representation can be constructed from 
the corresponding constrained integrable hierarchies. The KP hierarchy with self-consistent 
sources (KPHSCS) and the modified KP hierarchy with self-consistent sources (mKPHSCS) 
are constructed in this way and their interesting solutions are obtained via generalized Dar- 
boux transformation [26, 27]. In this sense, the soliton hierarchies with self-consistent sources 
may be viewed as integrable generalizations of the original soliton hierarchies. 

In [25], at the first time we studied the quasiclassical limit of the symmetry constraint 
of the Sato operator for the KP hierarchy in the framework of Sato theory which offers us a 
new way to study the reduction of the dispersionless hierarchy. Though the idea in [25] is 
different from that in [16], the well-known Zakharov reduction of the Sato function for the 
dispersionless KP (dKP) hierarchy can be also obtained as a result. The critical step in the 
procedure of this quasiclassical limit is to find the asymptotic forms of the potentials ) 
and rj(— ) appearing in the constraint of the Sato operator. Starting from the constrained 
Sato function for the dKP hierarchy, the constrained dKP hierarchy (cdKPH) and the dKP 
hierarchy with self-consistent sources (dKPHSCS) are therefore constructed. In contrast to 
the dKP case, the reduction of the dispersionless mKP hierarchy has been less studied. In 
this paper, we will further develop the idea in [25] to investigate the quasiclassical limit of 
the squared eigenfunction symmetry constraint of the Sato operator for the mKP hierarchy 
[31]. As will be shown below, this quasiclassical limit will lead to a reduction of the Sato 
function for the dispersionless mKP hierarchy. In this case, it is also important for us 
to find the asymptotic form of c[i(^) and rj(-) appearing in the constraint of the Sato 
operator for the mKP hierarchy. The constrained dmKP hierarchy (cdmKPH) will soon 
be obtained and its commutability can be proved. We can see that the dmKP hierarchy 
of (2 + l)-dimensions is decomposed into two commutative (1 + 1) -dimensional hierarchies 
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of hydrodynamical type. The dmKP hierarchy with self-consistent sources (dmKPHSCS) 
and its associated conservation equations will also be obtained according to our approach. 
Utilizing the conservation equations, we can develop the hodograph reduction method [6, 7] 
to solve the dmKPHSCS. In this sense, the dmKPHSCS may be regarded as an integrable 
generalization of the dmKP hierarchy. 

The paper will be organized as follows. In section 2, we briefly review some defini- 
tions and results about the mKPHSCS. In section 3, we take the quasiclassical limit to the 
squared eigenfunction symmetry constraint of the Sato operator for the mKP hierarchy. The 
cdmKPH and the dmKPHSCS together with its associated conservation equations are also 
obtained. Section 4 is devoted to presenting some solutions for the dmKPHSCS obtained by 
means of the hodograph reduction methodand. 

2 The squared eigenfunction symmetry constraint of 
the mKP hierarchy and the mKP hierarchy with self- 
consistent sources 

We first review some results about the mKP hierarchy with self-consistent sources in the 
framework of Sato theory [25-29]. Given a pseudo-differential operator (PDO) of the form 
[30] 

L = d + v + v 1 0- 1 + v 2 d~ 2 + (2.1) 
where Vi = fj(t), the mKP hierarchy is defined as 

d tk L = [Q k ,L], keN, (2.2) 

where Qk = {L k )>\ means the part of order> 1 of L k . The Lax euqation (2.2) is equivalent 
to the compatibility of the following linear equations 

L^mKP = tymKP, (2.3a) 
d tk TpmKP = Qk^mKP- (2.3b) 

i>mKP is often called the eigenfunction of the mKP hierarchy and also satisfies 

d\ipmKP = M mK pip mKP , (2.4) 

where M mKP is the Orlov operator of the mKP hierarchy. The adjoint eigenfunction i)* mKP 
satisfies 

L*lp* mK P = MmKP, ( 2 - 5a ) 
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dkVmKP = -QlVmKP, ( 2 -5b) 

btfnKP = -M* mKP r mKP - (2.5c) 
Making a constraint of the PDO L (2.1) as [31, 27] 

N 

L n = (L n )>! + ^ q^d-'r^d, n G N, (2.6) 
i=i 

where and rj(t) satisfying 

qi,t k = Qk(qi), r i>tk = -(dQ k d- 1 Y(r i ) = -{d^Qld)^), i = 1, AT, fc e N, (2.7) 

and Q fc = [(L n )«]>i = [((L n )>i + ^ili Qii^d^r^d) »]>i, we will get the (n-) constrained 
mKP hierarchy (cmKPH) as 

(L B ) tfc = [Q*,L B ], (2.8a) 

qi,t k = Qk(qi), (2.8b) 

r^-OT^Xr,), i = l,...,iV. (2.8c) 

(2.6) with (2.7) is called the squared eigenfunction symmetry constraint for the mKP hier- 
archy [31] and is called the integrated adjoint eigenfuntion, i.e., r,^ x satisfies (2.5b). 
If we add the term (Qk)t n to the right side of equation (2.8a) and requiring k < n, the mKP 
hierarchy with self-consistent sources (mKPHSCS) will be obtained as [27] 

(Q k ) tn -(L n ) tk + [Q k ,L n ] = 0, (2.9a) 

qi,t k = Qk(qi), (2.9b) 

^ = -(9-^)^), i = 1, N. (2.9c) 

As the case of KPHSCS [26] and many other cases in (l+l)-dimensions (see [22, 23] and 
the references therein), the cmKPH (2.8) may be considered as the stationary one of the 
mKPHSCS (2.9) if is viewed as the evolution variable. 

When k = 2, n = 3 in (2.9), we will get the mKP equation with self-consistent sources 
(mKPESCS) (with y = t 2 , t = t 3 ,v = v )[27] 

N 

13 3 3 

v t - -^Vxxx - -D'^Vyy) - -D~ l {v y )v x + -v 2 v x + ^2(qiri) x = 0, (2.10a) 

i=i 

ft, y = qi,xx + 2vq iiX , (2.10b) 
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n, y = -r itXX + 2vr itX , i = 1, • • • , N, (2.10c) 

where D^ 1 = f x is the integral operator. Under (2.10b) and (2.10c), (2.10a) will be obtained 
by the compatibility of the following auxiliary linear equations (Lax pair) 



i/> y = i/>xx + %vij> x , (2.11a) 

3 _ N 

•tyt = ipxxx + 3vip xx + -[v x + v 2 + D x 1 (v y )]ip x + S ^2q i D~ 1 (r i i) x ). (2.11b) 



i=i 



3 The quasiclassical limit 

In this section, we will consider the quasiclassical limit of (2.6) with (2.7) which will give 
rise to the dispersionless counterpart of (2.8) and (2.9), i.e. the constrained dispersionless 
mKP hierarchy (cdmKPH) and dispersionless mKP hierarchy with self-consistent sources 
(dmKPHSCS). 

Taking T n = et n and thinking of Vi(—) = Vi(T) + 0(e) as e — > 0, L in (2.1) changes into 

OO rji OO 

L^ed + ^Vii-^ed)- 1 = ed + Y,(Vi{T) + 0(e))(ed)-\ d = d x , X = ex. (3.1) 

i=0 i=0 

The constraint (2.6) now changes into 

N 

L^ = Q m + Y t q i (-)(ed)- 1 r i (-)ed, Q m = (L»)> 1 , (3.2) 



=i e e 



where <fc(j) and r^) satisfy 

tM-)]T k =QMj)), ^(-^^-[(ea)-^:,^)]^^-)), Q tk = [(L^]> 1 , i=l,...,JV. 

(3.3) 

It is easy to prove that 

OO 

C = ^(L e )=p + J2^(T)p- 1 - (3.4) 

i=0 

is a solution of the dmKP hierarchy, i.e., satisfies 

d Tn C = {Q n ,C}, (3.5) 

where a e denotes the principal symbol [9], the Poisson bracket is defined as 

r A , N , ^ dAdB dAdB 

Mfe ,), flta ,) } = __-__, (3.6) 
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and Q n = (C n )>i now refers to powers of p with order > 1. 

The dmKP hierarchy can be also written in the zero-curvature form 

^r-^f + {Q l ,Q s } = 0, I, seN. (3.7) 

When / = 2, s = 3, we will get the dmKP equation (Y = T 2 , T = T 3 and V = V ) 

2V T - *-D x \V YY ) - W X D-\V Y ) + 3V 2 V X = 0. (3.8) 



[32] and ipmKp(~) has the following WKB asymptotic expansion as e — > 0, 



Similarly like the dKP case [9], from (2.3), (2.4) and (2.5), it can be proved that ip m Kp(^) 



ipmKp(^, A) = exp[^S mKP (T, A)+0(1)], ^* mKP (j, A) = exp[~S mKP (T, A)+0(1)], e -> 0. 

(3.9) 

From (2.3b) and (3.9), we can obtain a hierarchy of conservation equations for the momentum 
function p = 8 ^ p , 

the compatibility of which, i.e., J^J^, = qj> dTi i m plies the dmKP hierarchy (3.7). 
Since qi{t) and ri(t) x evolve in the same rule as i) m KP and r ^* mKP respectively, we regard 

qi(J-) = ip mK p(^-, A = Aj) ~ exp{ SmKP< ^' — + a^i + 0(e)], (3.11a) 

e[r<(|)U = Ckp(^ A = A,) ~ exp[- SmKP ^ ^ + a t2 + 0(e)], e - 0, i = l,...,N. 

(3.11b) 

By a tedious computation, we will find that when e — > 0, 

)(e0) = + (^P(r, A,)x + 0(e))(69)- 1 + • • • 

+ ((S mi , P (T, A,)x) n + 0(e))(ea)- n + •••]. 

Taking the principal symbol of both sides of (3.2), we obtain the constraint for C as 

N i 

C n = Qn-T^ 77FT ^^ + S m Kp(T,X l ) x p- 1 + --- + (S mKP (T,\ l ) x ) n p" n + ---} 

N 

= &-£(- + —)' (3.12) 



i=i 

where Q n = (£ n )>i and 



a . = e a 11+ a l2) v . = SmKP (T, X t ) x . (3.13) 
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From (3.3), (3.11) and (3.13), a direct tedious computation shows the evolution for a« and 
Pi, i.e., 

d 

ai, Tk = [ai(^Qk(p))\ P = Pl }x, (3.14a) 
Pi, Tk = [Qk(p)\ P = Pi ]x, keW = l,...,N, (3.14b) 

where 

N 

Q k = [(£»)*]>! = {[Q n - + r^r)l-}>i- (3-15) 

<=1 ^ " 

It is easy to prove that with (3.14), the constraint (3.12) is compatible with the dmKP 
hierarchy (3.5). So the quasiclassical limit of the constrained mKP hierarchy (2.8) gives rise 
to the following constrained dispersionless mKP hierarchy (cdmKPH) 

(C n ) Tk = {Q k ,C n }, (3.16a) 
d 

ai > T k = [ a i(-Q^Qk(p))\ P = Pl }x, (3.16b) 

Pi,r k = [Qk(p)\ P = Pi ]x, i = 1,~,N, (3.16c) 
where C n and Qk are given by (3.12) and (3.15) respectively. 

Theorem 3.1 The equations with different k in (3.16) are commutative. 
Proof.From (3.16a), we can prove that for V fci^G N, the following identity holds, 

(Q kl ) Tk2 - (Q k2 ) Tkl + {Q kl ,Q k2 } = o, (3.17) 

k ■ 

where Q ki = [(L n )~^]>i, i = 1,2. A direct computation shows that 

{Pi)T kl T k2 = {[Q k i(p)\p= Pi \x}T k2 = [(Qk l )T k2 \p=p i + (Qk l )p\p=p i ((Qk2)x\p=p i + (Qk 2 )p\p=PiPi,x))}x- 

So, (pi) Tkl T k2 = (Pi)r k2 T kl is equivalent to 

{Q k ±)T k2 \p= Pi + (Q kl ) P (Qk 2 )x\ P = Pi = {Q k2 )T kl \ P = Pi + (Q k2 ) P {Q kl )x\ P = Pi , 

which holds from (3.17). 

Analogously, we can prove (a^^ = (a;)^^. The proof for (£ n ) Tfcl T fc2 = (£ n k 2 T fcl is 
similar to that for the mKP hierarchy. 
This completes the proof. 

In fact, the commutativity for the cmKPH (3.16) is a result of that for the mKP hierarchy 
since the constraint (3.12) and (3.14) is compatible with the mKP hierarchy. 
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As will be shown below, the (2 + l)-dimensional dispersionless mKP hierarchy (3.7) is de- 
composed into two commutative (1 + 1) -dimensional hierarchies of hydrodynamic type of 
(3.16) with k = s and k = I respectively. For example, when n — 1, (3.12) becomes 

N 

£ = p-J2(- + —)> ( 3 - 18 ) 

~i Pi p-Pi 

and we immediately get the following two systems of hydrodynamic type when k = 2 and 
A; = 3 (Y = T 2 , T = T 3 ), 

a iX = 2[oi(pi - Si)]*-, (3.19a) 
Pi , Y = \p1 - 2S lPi ] x , (3.19b) 

and 

a t , T = [a t (3p* - 6S lPi + 3Sf - 3S 2 )] X , (3.20a) 

Pi,T = [p! - 3S lP i + (3Sf - 3S 2 ) Pl ] x , (3.20b) 

where Si = ~ an< ^ ^ = YliLi a i- ^ is easy to verify that common solution of the 

systems (3.19) and (3.20) generates a solution V = — YliLi ~ °f the dmKP equation (3.8). 
Generally, if C n = Q n — Yli=i{^ + satisfies (3.16) with k = s and k = I simultaneously, 
then Qi = [(C n )«] > 1 with i — l,s will satisfy (3.7). 

When adding the term {Qk)T n to the right hand side of (3.16a) and requiring k < n, or 
taking the quasiclassical limit of (2.9) directly, we will obtain the dmKP hierarchy with 
self-consistent sources (dmKPHSCS) as 

(Q k ) Tn -(C n )T k + {Q k ,C n } = 0, (3.21a) 
d 

= [ a i(-Q^Qk(p))\ P = Pi }x, (3.21b) 

Pi,r k = [Qk(p)\ P = Pi ]x, i — 1, N, (3.21c) 

where C n and are given by (3.12) and (3.15) respectively. If "T„" is viewed as the 
evolution variable, the cdmKPH (3.16) may be regarded as the stationary system of the 
dmKPHSCS (3.21) like the dispersionfull case. It is not difficult to prove that under (3.21b) 
and (3.21c), the equation (3.21a) will be obtained by the compatibility of the following 
conservation equations 

p n = [Qk(p)]x, (3.22a) 

N 

PT n = [C n (p)]x = [Q n (p) ~ + -3-)U- (3.22b) 

i=l Pi P Pi 
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For example, when k = 2, n = 3, 

N 



C 3 = p 3 + 3l> 2 + (3V* + 3V 1 )p - + -^-), 

^ P; p-p* 



i=i 



(3.21) becomes the dmKP equation with self-consistent sources (dmKPESCS) (Y = T2, 
T = T 3 , V = V and V 1 is eliminated by V liX = \V^ Y ~ U V o)x) 

3 N 
2V T - -D- X \V YY ) - 3V X D X \V Y ) + 3V 2 V X - 2 = 0, (3.23a) 

1 i=i ^ 

Oi.y = + V)] x , (3.23b) 

p l ,y = (p l 2 + 2\/p J ) x , i = l,...,iV. (3.23c) 

Under (3.23b) and (3.23c), (3.23a) will be obtained by the compatibility of the following 
conservation equations 

p Y = 2 PPx + 2V xP + 2Vp x , (3.24a) 
PT = 3V xP * + (3V 2 + V 1 ) x p+(3p 2 + 6Vp + 3V 2 + V 1 ) Px -f2i(^) X + ^ 

j =1 Pi P Pi \P Pi) 

(3.24b) 

where V ljX = \V Y - \{V 2 ) X . 

As will be shown in the next section, utilizing the conservation equations (3.22), (3.21) 
will be solved by the hodograph reduction method. In this sense, the dmKPHSCS (3.21) 
can be regarded as an integrable generalization of the dmKPH (3.7). 



4 Hodograph reduction solutions 

In this section we will present some solutions of the dmKPHSCS (3.21) obtained by means 
of the hodograph reduction method [6, 25]. Provided the conservation equations (3.22), one 
can consider the M-reductions of (3.22) so that the momentum function p, the auxiliary 
potentials V^, % > 1 and a^, p iy i = 1, iV depend only on a set of functions W —(W 1 ,...,Wm) 
with Wi = Vo = V and (Wi,...,Wm) satisfy commuting flows 

dW dW 

wr MW) w n> - 2 

where the M x M matrices A n are functions of (W±, Wm) only. In the following, for 
example, we list some results for the dmKPESCS (3.23) in the cases of M — 1 and M — 2. 
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1. M = 1 

In this case, we will get the following hodograph equation 

,4c- 4 _ r 3 / 4c-4 s0 _ 3c 



X + C-^V)Y + ) V 2 + — ^- V 2 + 6 + E *)V^]T = F(V), (4.2) 

^ ' i=i 



where c,di are constants and F(V) is an arbitrary function of V. 
(l)c = 3, ^§ = 1, F(F) = 0. Then 



„ "[(Ef=i *)T + m ± ^[(Eli + 81? - 210(6T + X) 

V= , (4.3) 

105T ' v ; 

and V, Pi = 2V, a { = diV 3 , % = 1, • • • , N is an explicit solution for the dmKPESCS (3.23). 
(2)c = |, = 2, Choosing F(V) = 0. Then 



_ ± - (Ef =1 * + 81)(ftT + X)T 

(Eii* + 81)T 

and V, Pi = 3V, <n = d^V 4 , i = 1, • • • , TV is another solution for the dmKPESCS (3.23). 
Remark: When J2iLi di = 0, (4.3)) and (4.4)) will degenerate to the solution for the dmKP 
equation (3.8). 
2. M = 2 

In this case we can get the following solution of (3.23) as 

-Y ± J-Y* + 2(3T + 2) [(Ef =1 Ci)T + X] 
V = " 2(3T^) ' (45a) 

a . = _ Ci y(y + _J_), (4.5b) 

Pl = ~WT2~ v ' i = 1 >--> N - ( 45c ) 

with Cj are constants. 

Remark: (4.5a) degenerates to the solution of the dmKP equation (3.8) when Y^l=i c « = 0- 
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